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I.  INTRODUCTION 


The  analysis  of  flow-through  rocket  motor  exhaust  nozzles  has  under¬ 
gone  continuous  development  for  many  years,  since  the  optional  design  of 
these  nozzles  is  dependent  on  accurate  knowledge  of  the  flow  behavior  and  is 

important  to  the  attainment  of  high  thrust  efficiencies  for  launch  vehicles. 

1  2 

The  classic  analytical  solution  technique  based  on  the  series  expansion  ’  has 
limited  application,  as  it  requires  the  nozzle  entrance  to  be  suitably  shaped. 
During  the  past  decade  the  use  of  computers  for  the  solution  of  nozzle  flow 
fields  has  been  very  popular  among  research  engineers,  mainly  because  the 
modern  high-performance  propulsion  system,  for  the  sake  of  length  and  weight 
reduction,  usually  possesses  a  nozzle  contour  with  a  small  throat  radius  of 
curvature,  a  very  steep  wall  gradient  in  the  entrance  region,  or  a  submerged 
configuration,  and  the  numerical  technique  is  well-suited  for  application  to 
different  nozzle  geometric  configurations.  For  gas-only  one-phase  nozzle 
flows,  various  numerical  methods  used  in  the  past  were  reviewed  in  Ref.  9. 

^Hopkins,  D.  E.  and  D.  E.  Hill,  "Effect  of  Small  Radius  of  Curvature  on  Tran¬ 
sonic  Flow  in  Axisymmetric  Nozzles,  "  A  LA  A  J.,  4(8),  Aug.  1966,  p.  1337. 

2 

Kliegel,  J.R.  and  V.  Quan,  "Convergent-Divergent  Nozzle  Flows,  11  A  LA  A  J., 
Sept.  1968,  p.  1728. 

3 

Prozan,  R.  J.,  reported  in  "Numerical  Solution  of  the  Flowfield  in  the  Throat 
Region  of  a  Nozzle,  "  by  L.M.  Saunders,  BSVD-P-66TN-001  (NASA  CR82601), 
Aug.  1966,  Brown  Engineering  Co.,  Huntsville,  Ala. 

4 

Migdal,  D.,  K.  Klein,  and  G.  Moretti,  "Time- Dependent  Calculations  for 
Transonic  Nozzle  Flow,"  ALA  A  J.,  7^(1),  Feb.  1969,  p.  372. 

^Wehofer,  S.  and  W.  C.  Moger,  "Transonic  Flow  in  Conical  Convergent  and 
Convergent -Divergent  Nozzles  with  Nonuniform  Inlet  Conditions,  "  AIAA 
Paper  No.  70-635. 

^  Laval,  P.,  "Time -Dependent  Calculation  Method  for  Transonic  Nozzle  Flows," 
Lecture  Notes  in  Physics,  Jan.  1971,  p.  187. 

7 

Serra,  R.A.,  "Determination  of  Internal  Gas  Flows  by  a  Transient  Numerical 
Technique.  "  AIAA  J.,  10(5),  May  1972. 

Q 

Cline,  M.C.,  "Computation  of  Steady  Nozzle  Flow  by  a  Time -Dependent 
Method,  "  AIAA  J.,  12(4),  Apr.  1974,  p.  419. 

Q 

Brown,  E.  F.  and  G.  L.  Hamilton,  "A  Survey  of  Methods  for  Exhaust- Nozzle 
Flow  Analysis,  "  AIAA  Paper  No.  60,  1975. 


For  the  solid  rocket  motor,  one  of  the  prime  causes  of  performance 
loss  and  surface  damage  is  the  presence  of  condensed  metallic  oxide  particles 
of  the  combustion  products  in  the  flow  field.  The  thermal  and  velocity  lag 
associated  with  the  particles  often  results  in  decreased  nozzle  efficiency  and 
degradation  of  the  motor's  effectiveness  in  converting  from  thermal  to 
kinetic  energy.  Hence,  knowledge  of  the  role  played  by  the  nongaseous  com¬ 
bustion  products  in  the  rapid  expansion  through  the  throat  region  and  the 
qualitative  estimation  of  this  influence  are  essential  in  the  design  of  a  thrust 
nozzle.  A  comprehensive  review  of  investigations  involving  gas-particle 
nozzle  flow  fields  before  1962  is  presented  in  Ref.  10.  More  recent  studies 
include  the  numerical  iterative  relaxation  technique  of  Ref.  11  and  an  uncoupled 
flow  model  described  in  Ref.  12,  While  the  analysis  used  in  these  studies  is 
helpful  in  explaining  some  of  the  physical  processes  involved  in  the  gas- 
particle  flows  in  the  transonic  region,  both  suffer  from  the  same  weakness; 
i.e.,  the  assumption  that  the  gas-phase  streamline  coordinates  are  unaffected 
by  the  presence  of  particles.  This  assumption  is  particularly  inappropriate 
for  a  nozzle  with  a  very  small  throat  radius  of  curvature  or  very  steep  wall 
gradient,  since  the  presence  of  particles  can  alter  the  gas  flow  behavior.  The 
constant  fractional  l.'.g  and  the  linear  velocity  profile  assumptions  used  in 
Ref.  13  are  not  justified  a  priori.  The  results  obtained  or  refined  from  a 


^°Hoglund,  R.  F.,  "Recent  Advances  in  Gas-Particle  Nozzle  Flows,  "  ARS 
J ournal,  May  1962,  p.  662. 

^Regan,  J.  F.,  H.  D.  Thompson,  and  R.F.  Hoglund,  "Two-Dimensional 
Analysis  of  Transonic  Gas-Particle  Flows  in  Axisymmetric  Nozzles,  11 
J.  Spacecraft,  8(4),  Apr.  1971,  p.  346. 

12 

Jacques,  L.  J.  and  J.A.M.  Seguin,  "Two-Dimensional  Transonic  Two-Phase 
Flow  in  Axisymmetric  Nozzles,  "  ALAA  Paper  No.  74-1088,  Oct.  1974. 

13 

Kliegel,  J.  R.  and  G.  R.  Nickerson,  "Axisymmetric  Two- Phase  Perfect 
Gas  Performance  Program,  "  Report  02874-6006-R000,  Vols.  land  II, 

Apr.  1967,  TRW  Systems  Group,  Redondo  Beach,  Ca.  90278. 
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similar  analysis  for  the  transonic  region^  are  highly  uncertain,  although 
they  are  the  most  widely  used  method  in  the  propulsion  industry.  The  one¬ 
dimensional  analysis  shown  in  Refs.  15  and  16,  found  useful  in  some  areas, 
is  not  applicable  to  the  study  of  a  nozzle  with  a  steep  entrance. 

In  this  report,  the  time-dependent  method  is  applied  to  the  solution  of 

gas -only  one -phase  flow  and  fully  coupled  gas -particle  two-phase  flow  inside 

nozzles  of  arbitrary  geometry.  The  finite  difference  scheme  and  the  inlet 

boundary  conditions  incorporated  into  the  flow-field  program  are  shown  to 

yield  good  resolution  of  the  entire  subsonic-transonic  -  supersonic  flow  region. 

Moreover,  to  eliminate  the  computational  difficulty  associated  with  a  nozzle 

with  very  steep  wall  or  of  a  submerged  configuration,  the  Boundary- Fitted  - 

1 7 

Coordinates  (BFC)  system  is  adopted  for  generating  a  natural  grid.  Appli¬ 
cation  of  the  BFC  system  to  the  nozzle  flow  study  has  greatly  enhanced  the 
capability  of  the  flow-field  program  to  solve  problems  which  hitherto  have 
not  been  extensively  studied.  The  emphasis  of  the  study  has  been  placed  on 
one-  and  two-phase  flow  in  the  transonic  region.  Various  particle  sizes  and 
particle  mass  fractions  have  been  investigated  in  the  two-phase  flow.  The 
salient  features  associated  with  the  two-phase  nozzle  flow  compared  with 
those  of  the  one-phase  flow  are  illustrated  through  calculations  for  a  JPL 
nozzle  configuration,  for  the  Titan  III  solid  rocket  motor  nozzle,  and  for  the 
submerged  nozzle  configuration  utilized  in  the  IUS  solid  rocket  motor. 


14 

Coats,  D.  E.,  et  al.,  "A  Computer  Program  for  the  Prediction  of  Solid 
Propellant  Rocket  Motor  Performance,  11  Vols.  I,  II,  and  III,  AFRPL- 
TR-75-36.  July  1975. 

1 5 

Soo,  S.  L.,  "Gas  Dynamic  Processes  Involving  Suspended  Solids,  "  A.  I.  Ch.  E. 
Journal,  7(3),  Sept.  1961,  p.  384. 

1  6 

Hultberg,  J.A.  and  S.  L.  Soo,  "Flow  of  a  Gas-Solid  Suspension  Through  a 
Nozzle,  "  ALAA  Paper  No.  65-6,  Jan.  1965. 

17 

Thompson,  J.F.,  F.  C.  Thames,  and  C.W.  Martin,  "Boundary- Fitted 
Curvilinear  Coordinates  Systems  for  Solution  of  Partial  Differential  Equa¬ 
tions  on  Fields  Containing  Any  Number  of  Arbitrary  Two-Dimensional 
Bodies,  "  NASA  CR  2729,  July  1977. 
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II.  GOVERNING  EQUATIONS 


The  usual  assumptions  are  employed  below  to  derive  the  governing 
equations  of  a  gas-particle  two-phase  flow. 

a.  Mass  conservation  is  applied  to  both  mixture  and  individual 
phases. 

b.  The  mixture  flow  is  adiabatic,  i.  e.,  the  total  energy  of  the  mix¬ 
ture  is  constant. 

c.  Gas  phase  is  inviscid  except  for  its  interaction  with  the  metallized 
part  ^les,  where  the  momentum  exchange  is  considered  for  a 
viscous  gas  flow  over  spherical  condensed  particles. 

d.  Energy  exchange  between  the  gas  and  particles  occurs  through 
both  the  convective  and  radiative  heat  transfer. 

e.  The  particles  do  not  interact  with  each  other,  and  the  collision, 
condensation,  and  decomposition  of  the  particles  do  not  take  place. 

f.  The  gas  is  a  perfect  gas  and  is  chemically  frozen. 

g.  Volume  occupied  by  the  solid  particle  phase  is  negligible. 

Based  on  the  above  assumptions  and  normalized  by  the  gas-phase  stag¬ 
nation  state  corresponding  to  the  condition  at  the  inlet  plane,  the  governing 
equations  written  in  divergence  form  for  an  unsteady- state  two-phase  flow 
take  the  following  form: 


dE  dF  dG 
■§F  +  ^ 


+  H 


0 


(1) 
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where 


N  =  1  one -phase  gas -only 
N  =  2  two-phase  flow 

6  =  0  two -dimension 
6  =  1  axisymmetry 


The  nondimensional  parameters  used  here  are  gas -phase  stagnation  pressure 
Ptj,  stagnation  density  Ptj,  stagnation  energy  per  unit  volume  efcj  = 

Ptj/(T-1),  maximum  speed  and  stagnation  temperature  evaluated 

at  the  inlet  plane,  so  that 


P  =  P/Ptl 


r  *  (y-i)/2y  , 


t  =  v  ,t/L 

maxi 


P  -  P/Pn 

U  a  u/V 


maxi  ’ 


v  =  v/V 


maxi  ’ 


e  =  e/e 


tl 


pi 1 


u.  *  u./V  .  , 
j  j  maxi  ’ 


v.  =  v./V  ,  , 
j  j  maxi  ’ 


h.  =  h,/e.  , 
J  J  tl 


x  =  x/L. 

r  =  7/L 


for  L  a  reference  length  scale 
(e,  g.,  unit  foot) 


where  p,  u,  v,  p,  and  e  are  the  dimensioned  density,  horizontal  x-component 

velocity,  vertical  r-component  velocity,  pressure,  and  energy  per  unit 

volume,  respectively,  for  the  gas  phase;  and  p.,  u.,  v.,  and  h.  are  the  dimen- 

J  J  J  J 

sioned  density,  x-component  velocity,  r-component  velocity,  and  energy  per 
unit  volume,  respectively,  for  the  particle  phase.  There  are  also 


Friction  term: 


a  •  i  Vj  l 

>  2  m.72  V  , 

J  J  maxi 
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Energy  exchange  term: 


whe  re 


with 


t 

T 


T. 

J 

m. 

J 


r . 
J 


c 


P 


Cj 


a 


€  . 
J 

€ 


B. 


2r[q. 


-  *c<Ty 


T)  -  gr(€.T4  -  €  T4)] 


g  =  N  ,/6  f.P  , 
6  c  u j  j  r 


=  <7  r  .T J? / 3c  H  £. 
J  tl  p'g  j 


q.  •  zlq.  a  u.(u-u.)  +  v.(v-v.) 

J  J  J  3  J  J 

T.  =  T./T..  =  [h jyp.  -  (u2  +v2)]/w 
J  J  tl  J  3  J  J 

t  =  t/tu  =  p/p,  o>  =  cvr 


=  dimensioned  time 
s  dimensioned  gas  temperature 
=  gas  viscosity 

=  dimensioned  particle  temperature 

*  particle  mass  density 
=  particle  radius 

a  gas  specific  heat  at  constant  pressure 
■  particle  heat  capacity 
=  gas  specific  heat  ratio 
a  Stefan-Boltzmann  constant 
«  particle  emissivity 

*  gas  emissivity 
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r  s  dimensioned  radial  coordinate 


P  «  gas  Prandtl  number 
x  s  dimensioned  axial  coordinate 
The  momentum  transfer  parameter  f.  is  defined  as 


CD/CD 


Stokes 


where  Cn  is  the  drag  coefficient  based  on  C.  B.  Henderson's  correlation  equa- 
18  ^ 

tion  for  spheres  in  continuum  and  rarefied  flows,  and  Cr.  s  Re./24  is 

Stokes  J 

the  Stokes  law  of  drag  coefficient  for  spheres  in  creeping  motion. 

The  heat  transfer  parameter,  particle  Nusselt  number,  is  taken  as 


N  .  =  2  +  0.459  R°:  55  P°*  33 
uj  ej  r 


The  particle 

|dq.|/V 
I  j 1  max 

as  follows: 


Reynolds  number  is  based  on  the  relative  speed  |  il  q^  |  = 

.  =  >/(u-u.)^  +  (v-v.)^  and  the  particle  radius  r.,  and  is  < 
1  J  J  r  J 

l^i  I7/ 


defined 


R  .  * 
ej 


maxi 


The  gas  viscosity  is  evaluated  from 


**  =  » 


8  '•tl  - 


where /ltj  is  the  gas  viscosity  at  the  stagnation  temperature  T^  corresponding 
to  the  inlet  condition,  and  A  is  an  input  constant. 


18 


Henderson,  C.B.,  "Drag  Coefficients  of  Spheres  in  Continuum  and  Rarefied 
Flows,  "  A 1A A  J.,  _1_4(6),  June  1976,  p.  707. 
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It  is  often  debated  which  form  of  the  equation  the  drag  coefficient 
and  the  particle  Nusselt  number  Nu^  should  take.  The  correlating  equations 
of  Ref.  18  provide  accurate  representations  of  sphere  drag  coefficients 
over  a  wide  range  of  flow  conditions.  The  simple  form  of  the  particle  Nusselt 
number  is  adopted  from  Ref.  19.  If  more  advanced  parametric  equations  are 
available,  they  can  be  incorporated  easily  into  the  present  analysis  without 
much  modification. 

The  computer  program  Dusty  Transonic  Internal  Flows  (DTIF)  has  been 
developed  so  that,  for  the  gas-only  one-phase  flow  (N=l),  all  the  particle 
phase  calculations  are  bypassed. 


19 

Carlson,  D.  J.  and  R.  F.  Hoglund,  "Particle  Drag  and  Heat  Transfer  in 
Rocket  Nozzles,  "  A  LA  A  J.,  2jll),  1964,  p.  1980. 
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III.  NUMERICAL  ASPECTS 


From  a  general  arbitrary  nozzle  configuration  in  the  physical  plane  x,  r 
the  transformation  to  a  grid  with  uniform  square  mesh  in  the  computational 
plane  can  be  accomplished  by  using  BFC;  this  requires  the  solution  of 

17 

two  elliptical  partial  differential  equations  with  Dirichlet  boundary  conditions. 

Figure  1  illustrates  the  transformation  relationship.  The  solution  utilizing 

the  successive  over- relaxation  (SOR)  method  for  generating  the  boundary- 

fitted  coordinates  is  carried  out  by  the  TOMCAT  program,  and  the  scale 

17 

factors  for  transformation  are  computed  in  the  FATCAT  program. 

Formally  applying  the  chain  rule  of  change  of  independent  variables  for 
Eq.  (1)  results  in  the  following  conservation  laws  in  the  plane: 


dE  dF  dG 

Ti  +"Se 


+  H  =  0 


(2) 


whe  re 


E  s  EJ2 

F  =  Fr.  -  Gx» 
I  G  =  Fr^.  ^ 

H  =  HJ2 


and  J2  »  x 


x^r^.  is  the  Jacobian  of  transformation.  For  a  particular 


nozzle  geometry  and  transformation,  the  Jacobian  and  the  partial  derivatives 

j'j 

are  computed  in  the  FATCAT  program  and  stored  on  disk  for  flow  field  study. 


V 

For  a  simple  region  such  as  the  nozzle  geometry  considered  herein,  the 
FATCAT  program  gives  the  values  of  scale  factors  at  the  corner  points 
twice  as  large  as  they  should  be.  This  has  been  corrected  for  the  nozzle 
application  in  this  study. 


PHYSICAL  PLANE  lx,r) 


Through  insertion  of  the  unsteady  term  in  the  governing  Eq.  (1),  the 

differential  equation  is  cast  into  hyperbolic  type;  therefore,  the  complication 

associated  with  the  mixed  flow  phenomenon  necessarily  existing  in  a  steady- 

20 

state  analysis  is  eliminated.  The  MacCormack  finite  difference  scheme 

21. 

which  has  been  applied  successfully  to  nozzle  flow  problems  is  adopted  here 
for  the  solution  of  the  partial  differential  Eq.  (2). 

For  one -phase  flow,  the  initial  condition  is  based  on  a  one-dimens  ional 
isentropic  analysis  with  the  flow  vector  set  to  the  local  inclination  angle  from 
linear  interpolation  between  the  lower  and  upper  wall  slopes  along  the  same 
grid  line  (constant  £).  The  converged  one-phase  results  serve,  then,  as  the 
initial  guess  for  the  gas -phase  data  in  the  two-phase  flow. 

For  the  particle -phase  arrays  of  initial  velocity  lags  X and  temperature 
ratios  Xy  are  chosen,  and  the  initial  condition  (guess)  is 

p.  =  p<t>l(\-<t>),  u.  =  uXv,  v.  =  vXv 

T.  =  T/Xx,  h.  *  yp.[wT.  +  (uj  + 

where  <t>  -  W./W  is  the  particle  mass  fraction  and  0)  =  C./C  is  the  ratio  of 

j  m  J  P 

particle  heat  capacity  to  gas  specific  heat  at  constant  pressure. 

A  unit  velocity  lag  and  temperature  ratio  as  an  initial  guess  of  particle 
phase  are  satisfactory  for  this  study. 

The  initial  guess  based  on  the  so-called  Equilibrium  Gas-Particle 
Mixture  (EGPM),  which  is  itself  a  physically  meaningless  term  with  its  isen¬ 
tropic  index  derived  from 

v  =  y_L±w&k£L 

'm  1  +YU)<f>{l  -<f>) 

is  not  appropriate  for  the  fully  coupled  two-phase  flow  study. 

20MacCormack,  R.W.,  "The  Effect  of  Viscosity  in  Hypervelocity  Impact 
Cratering,  "  AIAA  Paper  69-354,  May  1969. 
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Chang,  I-Shih,  "Three-Dimensional  Supersonic  Internal  Flows,  "  AIAA 
Paper  76-423,  July  1976. 
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The  exit  boundary  condition  is  based  on  a  simple  linear  extrapolation 
since  the  mixture  flow  is  assumed  to  be  supersonic  at  the  exit  plane,  and  the 
error  generated  from  the  extrapolation  is  not  expected  to  propagate  back 
and  affect  the  upstream  results. 

The  inlet  boundary  condition  for  horizontal  inflow  is  computed  from  a 
characteristics  formulation  similar  to  that  of  Ref.  7,  which  provides  fairly 
smooth  subsonic  flow  in  the  physical  domain.  For  radial  inflow,  e.g.,  the 
IUS  motor  studied  herein,  the  experimentally  evaluated  propellant  burning 
rate  and  chamber  pressure /temperature  data  supply  needed  information  at 
the  propellant  burning  surface,  and  a  linear  interpolation  for  smoothing  flow 
variables  at  the  grid  line  adjacent  to  the  propellant  burning  surface  is  required 
to  avoid  instability  in  the  inlet  region. 

For  a  nozzle  with  a  centerbody,  the  flow  variables  at  the  boundary  are 
obtained  from  linear  extrapolation  of  the  data  from  two  adjacent  interior 
points  and  then  modified  by  the  local  tangency  condition.  Without  the  center- 
body  in  the  axisymmetric  nozzle,  the  flow  variables  take  undetermined  form 
at  the  centerline.  The  standard  L'Hospital's  rule  with  the  symmetry  con¬ 
sideration  is  used  to  evaluate  all  the  physical  flow  variables  except  the  radial 
velocity,  which  is  zero  at  the  singular  centerline.  Also,  if  notice  is  taken 
that  all  the  conservative  variables  at  the  centerline  are  zero  except  =  -rp 
(i.  e.,  no  restriction  is  imposed  on  the  gas  pressure  change  at  the  centerline), 
a  smoothing  process  involving  a  linear  interpolation  for  resetting  the  gas- 
phase  radial  velocity  v  at  the  first  interior  point  above  the  centerline  is 
found  helpful  in  stabilizing  the  solution. 

It  is  important  also  to  retain  the  fourth-order  damping  terms  to  the 
second-order  MacCormack  finite  difference  scheme  for  the  unsteady  state 
application  in  order  to  eliminate  the  nonlinear  instability.  The  formulation 
adopted  here  is  similar  to  that  of  Ref.  22  with  a  damping  coefficient  equal  to 
0.01. 

22  “  ‘ 

Kutler,  P.,  L.  Sakell,  and  G.  Aiello,  "On  the  Shock-on-Shock  Interaction 
Problem,  "  AIAA  Paper  No.  74-524,  June  1974, 
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The  integration  step  size  is  governed  by  the  local  CFL  condition 
similar  to  that  used  in  Ref.  5  and  determined  by  the  following  expressions 


where 


and  a  = 


At  =  min 


Al 
q  +  a 


A  l  =  'J  (Ax)^  +  {Ar)^,  q  =  /7T7" 
)/2  is  the  dimensionless  local  sonic  speed. 
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IV.  NOZZLE  WITH  SMALL  THROAT  RADIUS  OF 
CURVATURE --JPL  NOZZLE 


A.  ONE-PHASE  FLOW 

23 

The  compressible  flow  inside  the  JPL  axisymmetric  nozzle  with 
45°  entrance  and  15°  exit  straight  wall  tangent  to  a  circular  throat  (with  ratio 
of  throat  radius  of  curvature  to  throat  height  =  0.625)  provides  a  classic 
comparison  for  the  present  nozzle  flow  study.  Figure  2  shows  the  physical 
grid  generated  from  the  boundary -fitted  coordinates  system  mentioned  above. 
The  computed  Mach  number  distributions  along  the  wall  and  along  the 
centerline  are  illustrated  in  Fig.  3;  the  test  data  are  also  shown  for  compari¬ 
son.  Figure  4  is  the  Mach  number  contour  plot,  and  Fig.  5  shows  the  Mach 
number  at  all  the  grid  points  in  the  physical  plane.  The  theoretical  gas-only 
one -phase  result  from  this  study  agrees  very  well  with  the  test  data  in  the 
entire  subsonic-transonic- super  sonic  flow  region.  This  can  be  attributed  to 
the  good  resolution  of  the  boundary  flow  variables  through  the  use  of 
boundary -aligned  grid  arrangement.  Note  the  smoothness  of  the  Mach  num¬ 
ber  distribution  in  the  subsonic  region  computed  by  the  present  method. 

The  recompression  waves  in  the  supersonic  region,  which  necessarily  occur 
due  to  over -expans ion  of  the  flow  near  the  wall  downstream  of  the  throat  with 
small  radius  of  curvature,  eventually  coalesce  into  a  shock  wave  near  the 

centerline  in  the  far-downstream  region.  This  flow  behavior  has  been  ob- 
7 

served  in  the  test  and  is  visible  from  Figs.  4  and  5. 

The  convergence  criterion  used  in  all  the  calculations  shown  herein 
requires  that  the  difference  in  Mach  number  is  at  least  0.01%  and  in  the  mass 
flow  rate  is  0.001%  at  the  throat  for  three  consecutive  time  integration  steps. 
For  the  JPL  nozzle  with  61  x  31  grid  points,  the  converged  solution  requires 
623  integration  steps  and  takes  6  min,  17  sec  execution  time  on  the  CDC  7600 
computer.  The  theoretical  results  agree  well  with  the  test  data  for  this 
23 

Cuffel,  R.  F.,  L.H.  Back,  and  P.  F.  Massier,  "Transonic  Flowfield  in  a 
Supersonic  Nozzle  with  Small  Throat  Radius  of  Curvature,  "  ALAA  J.,  7(7), 
July  1969,  p.  1364. 


Fig.  2.  BFC  Grid  for  JPL  Nozzle 


Mach  Number  Distribution  at  Wall  and  Centerline  fo 
JPL  Nozzle  (One-Phase  Flow) 
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nozzle  study,  thus  assuring  further  application  of  the  computer  code  to 
other  nozzle  configurations. 

B.  TWO-PHASE  FLOW 

The  same  flow  field  program  is  applied  to  the  fully-coupled  gas-particle 
two-phase  nozzle  flow  with  the  two-phase  index  N  in  Eq.  (1)  set  to  2.  For 
the  two-phase  flow  calculation,  the  following  data  are  adopted. 

Gas  Phase  (Air)  Particle  Phase 

C  =  0.255  Btu/lb  - °R  C.  =  0.33  Btu/lb  -°R 

P  m  j  m 

jT,  -  1.8xl0"5  lb  /ft-sec  m.  =  250  lb  /ft3 

tl  m  j  m 

y  =  1.4,  Pr  =  0.  74 
A  =  0.6 

The  chamber  condition  is  taken  to  be  T^  =  1000°R,  P^  =  150  psia. 

Also,  =  0.  1  and  €  =  0.  05  are  used  for  the  radiative  heat  exchange  between 
the  gas  and  spherical  particles. 

The  previously  computed  one-phase  flow  result  is  taken  as  the  initial 

guess  for  the  two-phase  flow.  Different  particle  sizes  and  particle  mass 

fractions  W./W  are  calculated.  Figure  6  shows  the  variation  of  the  com- 
J  m 

puted  throat  Mach  number  along  the  wall  and  along  the  centerline  at  each 

timewise  integration  step  (iteration)  for  various  particle  sizes  at  the  same 

W./W  =  30%.  At  300  iterations  the  two-phase  transonic  flow  region  is 
j  m  r 

essentially  established.  Further  integration  does  not  produce  appreciable 
change  in  the  throat  flow  field,  as  is  evidenced  from  the  continued  calcula¬ 
tion  with  r.  =  IfJ  to  600  integration  steps.  Figure  7  shows  the  computed  wall 
and  centerline  gas  phase  Mach  number  distribution  for  various  particle  sizes 
with  the  same  W./W  =  30%.  For  comparison,  the  previously  computed 
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Fig.  6.  JPL  Nozzle  Throat  Mach  Number  at  Every 

Integration  Step  (Two- Phase  Flow 

W./W  =  30%) 
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Fig.  7. 


J PL  Nozzle  Mach  Number  Distribution  at  Wall  and 
Centerline  (Two-Phase  Flow  W./W  =  30%) 


gas-only  one-phase  results  are  plotted  as  dashed  curves.  Some  of  the 
features  associated  with  the  fully-coupled  two-phase  nozzle  flow  in  compari¬ 
son  with  that  of  gas-only  one-phase  flow  are  found  in  Figs.  6  and  7.  While 
lower  gas  speed  is  observed  both  on  the  wall  and  centerline  in  the  two-phase 
flow  field  than  that  in  the  gas-only  one-phase,  the  small-sized  particles  act 
more  effectively  to  slow  down  the  gas-phase  expansion  than  that  of  large¬ 
sized  particles  for  the  same  particle  mass  fraction.  This  is  physically 
correct,  since  for  the  same  particle  mass  fraction  the  total  particle  surface 
area  effective  for  momentum  and  energy  exchange  between  gas  and  particles 
is  greater  in  a  two -phase  flow  field  involving  smaller  diameter  particles. 
Figure  8  gives  the  velocity  lag  q^/q  at  the  throat  and  at  the  exit  plane, 
respectively,  and  Fig.  9  shows  the  results  for  the  temperature  ratio  T/T^. 

The  temperature  ratio  defined  here  illustrates  the  thermal  lag  between  the 
particle  and  gas  and  is  less  sensitive  to  the  small  variation  in  the  computed 
numerics  than  the  commonly  used  thermal  lag  (1  -T.)/(l-T),  especially  in 
the  low  subsonic  region  when  both  gas  and  particle  temperatures  are  very 
close  to  the  reference  gas  stagnation  temperature  at  the  inlet  plane.  The 
large  particle  in  the  two-phase  flow  field  lags  much  more  behind  the  gas 
phase  flow  and  exhibits  a  wider  region  of  particle  free  zone  than  that  of  the 
small  particle.  Moreover,  in  comparison  with  the  present  fully  coupled 
two-phase  analysis,  the  constant  fractional  lag  assumption  used  in  Ref.  13 
is  justified  only  for  a  two- phase  flow  with  high  loading  ratio  of  very  small 
particles . 

The  effect  of  different  particle  sizes  in  the  two-phase  flow  can  be 
visualized  best  by  comparison  of  the  calculated  particle  density  contour  and 
particle-phase  velocity  vector  plot  for  the  small  (r^  =  1  fi)  and  the  large 
(r^  s  20 n)  particles  depicted  in  Figs.  10  and  11.  These  figures  show  that  for 
the  flow  with  lfj  radius  particles  a  sharp  change  in  particle  density  is  obtained 
near  the  upper  wall  downstream  of  the  throat,  and  the  particle  density 
drastically  decreases  to  a  small  value.  But  it  does  not  vanish  exactly  at  the 
wall.  With  the  large  (r^  =  20/J )  particle  flow,  however,  a  very  distinctive 


Fig.  9.  Gas-to- Particle  Temperature  Ratio 
(Two-Phase  Flow  W./W  =  30%) 
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particle-free  zone  appears  in  the  calculated  result  and  can  be  seen  from 
either  Fig.  10  or  Fig.  11.  Since  the  heavier  particles  apparently  cannot 
effectively  turn  around  the  throat  corner  with  small  throat  radius  of  curva¬ 
ture  and  evidently  tend  to  cluster  near  the  centerline,  there  are  essentially 
no  particles  present  near  the  wall  downstream  of  the  throat  to  slow  down  the 
gas  expansion.  This  explains  why,  in  Fig.  7,  the  difference  between  the 
Mach  numbers  in  one -phase  and  in  large  particle  two- phase  flows  near  the 
nozzle  lip  region  is  virtually  nil,  but  not  so  at  the  exit  centerline  region. 

The  variation  of  the  computed  throat  gas  phase  Mach  number  along  the 

wall  and  along  the  centerline  at  each  timewise  integration  step  for  a  different 

particle  mass  fraction  W./W  at  a  given  particle  radius  r.  *  ill  is  indicated 

J  m  j 

in  Fig.  12.  Figure  13  shows  the  corresponding  wall  and  centerline  gas  phase 
Mach  number  distribution.  As  before,  the  dashed  curves  are  the  results  for 
the  gas-only  one-phase  flow.  It  is  obvious  from  these  figures  that  a  reduc¬ 
tion  of  particle  mass  fraction  immediately  reduces  the  two-phase  "drag" 
effect.  A  high  particle  mass  fraction  (W ,/W^  Si  45%)  produces  an  entirely 
subsonic  flow  at  the  geometric  throat  for  the  nozzle  geometry  considered. 

It  is  then  possible  to  adjust  the  exit  Mach  number  from  supersonic  all  the 
way  down  to  subsonic  by  varying  the  particle  mass  fraction  and/or  the  particle 

size.  The  lip  shock  extending  to  the  exhaust  plume  field  which  occurs  in  the 

24  2  5 

over-  and  under-expanded  case  ’  could  also  be  weakened  or  eliminated 
through  the  particle  drag  effect. 

A  close  look  at  Figs.  6  and  12  also  reveals  the  fact  that  care  must  be 

taken  in  the  case  of  large  particle  mass  fraction  flows  (W./W  >  45%)  and 

_  J  ^ 

very  small-sized  particle  flows  (r.  <  0.  5/i)  in  starting  from  the  one-phase 
initial  guess.  A  sharp  drop  in  the  wall  Mach  number  at  or  near  downstream 

24 

Chow,  W.  L.  and  I-Shih  Chang,  "Mach  Reflection  Associated  with  Over- 
expanded  Nozzle  Free  Jet  Flows,  "  AIAA  J.,  13(6),  June  1975. 
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Chang,  I-Shih  and  W.  L.  Chow,  "Mach  Disc  from  Underexpanded  Axi- 
symmetric  Nozzle  Flow,  "  AIAA  J.,  12,  Aug.  1974,  p.  1079. 
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Fig.  1 1.  JPL  Nozzle  Mach  Number  Distribution  at  Wall  and 
Centerline  (Two- Phase  Flow  r.  =  Ifl) 


of  the  throat  can  be  expected,  especially  when  the  throat  radius  of  curvature 
is  small,  necessitating  the  use  of  small  time  marching  step  at  the  beginning 
of  the  two-phase  calculation.  Furthermore,  the  possibility  of  subsonic  flow 
occurring  at  the  exit  plane  requires  modification  of  the  supersonic  down¬ 
stream  boundary  condition  or  the  extention  of  the  exit  boundary  to  a  farther 
downstream  location. 

Figures  14  and  15  show  the  computed  velocity  lag  and  temperature 
ratio  for  various  particle  mass  fractions  at  a  fixed  r.  =  Ifl.  Figure  16  sum¬ 
marizes  several  Mach  number  contours,  M  =  0.  1,  0.  2,  0.  5,  1. 0,  and  1.  5, 

for  various  particle  sizes  at  a  fixed  ratio  W./W  =  30%,  and  Fig.  17  is  the 

J  m  6  _ 

similar  result  for  various  particle  mass  flow  ratios  at  a  fixed  radius  r.  =  1^. 

J 

As  before,  the  dashed  curves  are  the  results  for  a  gas-only  one-phase  flow. 
A  typical  two-phase  run  involving  300  integration  steps  takes  approximately 
7  min  execution  time  on  the  CDC  7600  computer. 
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V.  NOZZLE  WITH  VERY  STEEP  ENTRANCE- -TITAN  III  MOTOR 


A.  ONE-PHASE  FLOW 

A  severe  test  of  the  present  numerical  technique  is  the  solution  of  the 
compressible  flow  field  inside  a  nozzle  with  near  vertical  entrance  region 
shown  in  Fig.  18,  which  also  illustrates  the  grid  generated  from  the  boundary- 
fitted  coordinates  system.  The  nozzle  steep  entrance  contour  is  that  of  the 
Titan  III  solid  rocket  motor  nozzle.  The  specific  heat  for  the  combustion 
products  is  Y  -  1.  19.  The  computed  wall  and  centerline  Mach  number  distri¬ 
bution  is  given  in  Fig.  19.  The  same  nozzle  was  analyzed  with  Y  =  1.4,  and 
the  results  are  also  plotted  on  the  same  figure  which  serves  to  illustrate  the 
effect  of  different  y  on  the  Mach  number  distribution.  Figure  20  depicts  the 
computed  throat  Mach  number  at  every  integration  step.  The  converged 
one-phase  flow  solution  for  Y  -  1.19  requires  890  integration  steps  and  takes 
8  min,  58  sec  on  the  CDC  7600  computer.  Figure  21  is  the  Mach  number 
contour  plot,  and  Fig.  22  shows  the  Mach  number  at  all  the  grid  points.  The 
sonic  point  on  the  wall  has  been  observed  to  be  far  upstream  of  the  throat, 
which  indicates  that  higher  heating  rate  occurs  farther  upstream  of  the  throat 
than  that  expected  from  the  simple  one-dimensional  analysis.  This  partially 
explains  why,  in  past  full-scale  firings,  the  Titan  III  motor  nozzle  was 
ablated  much  more  in  a  region  far  upstream  of  the  throat  than  at  the  throat. 
Cold-flow  tests  were  recently  conducted  at  the  Chemical  Systems  Division 
of  United  Technology  Corporation  for  the  6°  canted  Titan  III  solid  rocket 

motor  using  nitrogen  with  Y  -  1.4;  the  wall  Mach  number  at  the  throat  was 
27 

measured  to  be  1.  52  which  agrees  fairly  well  with  the  computed  value  1. 6 
for  the  axisymmetric  nozzle. 


"’Private  communication  with  Chemical  Systems  Division  of  United  Tech¬ 
nology  Corporation  on  Titan  III  Engineering  Drawings,  Nov.  1978. 


Private  communication  with  R,  Dunlap,  Chemical  Systems  Division  of 
United  Technology  Corporation,  Nov.  1978. 


Fig.  18.  BFC  Grid  for  Steep  Entrance  Nozzle 
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Fig.  19.  Mach  Number  Distribution  at  Wall  and  Centerline  for 
Steep  Entrance  Nozzle 


Fig.  20.  Throat  Mach  Number  at  Every  Integration  Step  for 
Steep  Entrance  Nozzle 


Fig.  21.  Mach  Number  Contour  for  Steep  Entrance  Nozzle 


58.991 


-53.1626 
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Fig.  22.  Mach  Number  Pictorial  Plot  for  Steep  Entrance  Nozzle 
(One-Phase  Flow) 


mm 


Note  that  no  converged  solution  is  possible  with  the  conventional  grid 
with  vertical  axial  stations,  such  as  that  used  in  Refs.  4  through  8  and  in 
most  of  the  nozzle  studies  reported  thus  far,  for  the  nozzle  with  a  steep  wall 
slope  like  that  of  the  Titan  III  motor.  It  is  the  author's  experience  that  when 
the  nozzle  wall  slope  is  greater  than  «60°,  the  conventional  grid  with  vertical 
axial  station  cannot  handle  the  drastic  change  in  flow  properties  along  the 
steep  wall,  and  this  results  in  numerical  instability.  On  the  contrary,  no 
difficulty  is  encountered  in  the  calculation  with  the  grid  generated  from  the 
boundary -fitted  coordinates  system. 

B.  TWO-PHASE  FLOW 

The  two-phase  flow  data  for  the  Titan  III  motor  are  as  follows: 


Gas  Phase 

C  =  0.  64  Btu/lb  -°R 
p  m 

ZT  .  =  5.  97x10  3  lb  /ft-sec 
^tl  m 


Particle 

Phase 

c. 

=  0.33  Btu/lb  -°R 

J 

m 

m. 

a  200  lb 

/ft3 

J  m 


P  r  0. 45  r.  =  6/i 

r  ) 

A  =  0. 664,  y  =  1.  19  W./W  =  28.  8% 

J  m 

The  chamber  condition  is  T^j  *  5890°R,  =  600  psia. 

The  variation  of  the  computed  throat  gas -phase  Mach  number  along  the 
wall  and  at  the  centerline  at  every  integration  step  based  on  the  initial  guess 
from  the  previously  computed  one-phase  results  is  also  illustrated  in  Fig.  20. 
At  the  end  of  500  integration  steps  the  gas -phase  Mach  number  distribution  is 
shown  in  Fig.  19  for  comparison  with  that  of  one-phase  flow.  The  gas-phase 
Mach  number  contour  is  plotted  in  Fig.  21  and  the  corresponding  gas -phase 
Mach  number  distribution  throughout  the  flow  field  is  depicted  in  Fig.  23. 
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The  particle  density  contour  and  distribution  are  plotted  in  Figs.  24  and  25. 
The  gas-phase  pressure  field  has  not  been  altered  much  by  the  introduction 
of  particles  in  the  flow  field,  and  the  pressure  distribution  is  shown  in 
Fig.  26.  The  velocity  and  temperature  ratio  are  indicated  in  Fig.  27.  The 
calculation  of  500  integration  steps  for  the  two-phase  flow  takes  10  min, 

3  7  sec  execution  time  on  the  CDC  7600  computer. 

Particles  are  likely  to  impinge  on  the  steep  wall  in  the  entrance  region. 
The  present  analysis  does  not  incorporate  any  pertinent  particle  impinge¬ 
ment  model  for  calculating  the  erosion  caused  by  such  impingement.  Never¬ 
theless,  the  particle  density  contour  obtained  from  this  study  does  show 
regions  of  high  particle  concentration  which  may  affect  the  results  from 
boundary  layer  calculations  and  thereby  the  results  of  transient  in-depth 
thermal  analyses  for  the  prediction  of  nozzle  wall  temperature. 
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Fig.  24.  Particle  Density  Contour  for  Steep  Entrance  Nozzle 
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Fig.  25.  Particle  Density  Pictorial  Plot  for 
Steep  Entrance  Nozzle 
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26.  Pressure  Distribution  for  Steep  Entrance  Nozzle 
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Fig.  27.  Velocity  Lag  and  Temperature  Ratio  for 
Steep  Entrance  Nozzle 


VI.  SUBMERGED  NOZZLE- -IUS  SMALL  MOTOR 


A.  ONE-PHASE  FLOW 

It  has  long  been  recognized  that  the  solution  of  the  flow  field  inside 
rocket  motors  with  a  submerged  nozzle  configuration  constitutes  an  impor¬ 
tant  phase  of  the  flow-field  study.  Both  large  and  small  IUS  solid  rocket 
2  8 

motors  have  a  submerged  nozzle.  The  IUS  will  be  used  as  an  upper  stage 
for  both  the  Titan  III  and  the  Shuttle.  In  the  past,  due  to  the  difficulty  of 
analyzing  the  internal  flow  field  for  IUS- like  motors  with  complicated  geo¬ 
metry,  various  approximations  have  been  adopted,  and  the  accuracy  of  the 
pressure  and  heat  transfer  predictions  on  the  submerged  nozzle  surface  was 
therefore  uncertain.  Although  the  viscous  effect  would  probably  dominate 
some  part  of  the  submerged  flow  region  in  the  gas-only  one-phase  flow,  the 
inviscid  flow  solution  shown  here  constitutes  a  first  attempt  toward  a 
complete  viscous  flow  solution  in  future  studies. 

The  IUS  small  motor  interior  geometry  including  the  igniter,  submerged 
nozzle  block,  and  a  propellant  burning  surface  is  illustrated  in  Fig.  28,  where 
the  physical  region  for  computation  has  been  identified  by  heavy  solid  lines. 
The  grid  generated  from  the  bounda ry- fitted  coordinates  system  is  depicted 
in  Fig.  29.  The  flow  region  is  bounded  by:  (a)  the  motor  case,  (b)  the 
motor  centerline  of  symmetry,  (c)  the  igniter  surface,  (d)  the  supersonic 
exit  plane,  and  (e)  the  propellant  surface  with  radial  mass  inflow.  The  region 
depicted  in  Fig.  29  incorporates  the  entire  subsonic  flow  region  without 
introducing  a  fictitious  vertical  inlet  boundary.  The  blownup  figure  for  the 
submerged  and  throat  region  is  shown  in  Fig.  30. 

The  propellant  burning  rate  for  the  small  IUS  solid  rocket  motor  is 
found  to  be  0.  206  in/sec  at  T^  =  5985°R  and  P  ^  =  410  psia,  which  results  in 

2  g 

"Inertial  Upper  Stage  SRM-II  Baseline  Design  Review,  "  Chemical  Systems 
Division  of  United  Technology  Corporation,  Dec.  1978. 


74.22  in. 


Fig.  28.  IUS  Small  Motor  Interior  Configuration  and 
Computational  Region 


AXISYMMETRY  xlin.l  x  =  5in. 


Fig.  29.  BFC  Grid  for  Small  IUS  SRM  with 
Submerged  Nozzle  Block 
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Blown -up  BFC  Grid  in  the  Submerged  and 
Throat  Region  for  Small  IUS  SRM 


l 

t 


the  fixed  inlet  condition  at  the  propellant  burning  surface  v  =  11.53  ft/sec 

(M  =  0.  00322)  for  C  =  0. 45  Btu/lb  -°R  and  y  =  1.  19.  The  Mach  number 

p  m 

in  the  junction  region,  where  the  igniter  joins  the  zero  radius  centerline, 
was  computed  incorrectly  in  an  earlier  study.  In  this  report,  the  boundary 
points  with  the  radial  coordinate  smaller  than  the  radial  length  of  the  adjacent 
finite  difference  mesh  are  treated  as  the  centerline  points  of  zero  radius,  and 
L'Hospital's  rule  is  conveniently  applied  thereby  avoiding  numerical  error 
resulting  from  decoding  conservative  variables  divided  by  a  very  small 
number  (small  radial  coordinate). 

Unlike  the  previous  two  nozzle  flows,  a  stricter  convergence  criterion 
is  deemed  necessary  for  the  submerged  nozzle  calculation  and  requires  that 
the  difference  in  Mach  number  be  less  than  0.  001%  and  in  mass  flow  rate 
less  than  0.001%  at  the  throat  for  three  consecutive  integration  steps.  For 
the  61  x  31  grid  points  shown  in  Fig.  29,  the  converged  solution  requires 
4487  integration  steps  and  takes  38  min,  16  sec  execution  time  on  the  CDC 
7600  computer. 

Figure  31  indicates  the  computed  throat  Mach  number  at  every  inte¬ 
gration  step.  The  IUS  solid  rocket  motor  has  a  throat  with  a  large  radius  of 
curvature;  the  computed  Mach  number  at  the  throat  is  1.071  on  the  wall 
and  0.  947  at  the  centerline  which  are  close  to  a  uniform  one-dimensional 
flow.  Figure  32  depicts  the  Mach  number  distribution  along  the  boundary, 
while  the  computed  Mach  number  contour  is  plotted  in  Fig.  33.  Figure  34 
shows  the  velocity  vector  plot  in  the  submerged  and  throat  region. 

13.  TWO-PHASE  FLOW 

The  following  data  are  used  for  the  two-phase  flow  inside  the  small 
IUS  motor: 
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Fig.  31.  Throat  Mach  Number  at  Every  Integration  Step  for 
Small  IUS  SRM 


Fig.  32.  Mach  Number  Distribution  on  the  Boundary  for 
Small  IUS  SRM  Nozzle 
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Fig.  34,  Velocity  Vector  Plot  in  the  Submerged  and  Throat 
Region  for  Small  IUS  SRM  (One-Phase  Flow) 


. . . 


Gas  Phase 

Particle  Phase 

c 

p 

=  0.  45  Btu/lb  -°R 
m 

C. 

3 

=  0.  32  Btu/lb  -°R 
m 

*tl 

=  5.  6  74x10"  3  lb  /ft-sec 
m 

m . 

3 

=  200  lb  /ft3 
m 

p 

r 

=  0.269 

W./W  =  30% 

J  m 

A 

=  0.65 

r . 

J 

=  2.5// 

y  =1.19 


The  chamber  condition  is  T^  =  5985°R,  P  .  =  410  psia. 

The  throat  Mach  number  at  every  integration  step  is  shown  in  Fig.  31 
for  easy  comparison  with  the  one-phase  solution.  The  Mach  number  distri¬ 
bution  along  the  boundary  surface  is  indicated  in  Fig.  32.  The  gas-phase 
Mach  number  at  exit  plane  is  1.57  at  centerline  and  2.41  at  wall  for  the  two- 
phase  flow;  these  are  smaller  than  the  corresponding  Mach  numbers  of 
2.25  and  2.  57  found  in  the  gas-only  one-phase  flow,  due  to  the  presence  of 
solid  particles  in  the  flow  field.  This  implies  that  an  IUS  solid  rocket  motor 
nozzle  flow  field  and  heat  transfer  analysis  based  solely  on  a  gas -only  one - 
phase  study  will  be  in  error.  Figure  35  shows  the  gas-phase  and  Fig.  36  the 
particle  velocity  vector  plot.  A  distinctive  particle-free  zone  is  visible  from 
the  particle  velocity  vector  plot  of  Fig.  36.  Figure  37  is  the  computed  gas- 
phase  Mach  number  contour,  and  Fig.  38  is  the  particle  density  contour  in 
the  two-phase  flow.  The  velocity  lag  and  gas -to-particle  temperature  ratio  at 
throat  (where  r  =  2.  15  in.)  and  exit  plane  (where  r^  =  4.05  in.)  are  shown  in 
Fig.  3  9. 

Although  the  submerged  nozzle  configuration  is  complex  and  the  govern¬ 
ing  two-phase  partial  differential  equations  are  highly  nonlinear,  no  compu¬ 
tational  difficulty  has  been  encountered  during  the  course  of  this  study,  and  the 
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Particle  Velocity  Vector  Plot  in  the  Submerged  and 
Throat  Region  for  Small  IUS  SRM  (Two- Phase  Flow) 


i 


timewise  integration  has  been  carried  out  in  a  straightfoward  manner. 

For  the  submerged  small  IUS  solid  rocket  motor  nozzle,  the  two-phase 
flow  field  calculation  of  1000  integration  steps  takes  23  min,  2  sec  execution 
time  on  the  CDC  7600  computer.  All  the  two-phase  flow  features  mentioned 
previously  for  the  JPL  and  Titan  III  nozzle  are  equally  applicable  to  this 
submerged  IUS  solid  rocket  motor  nozzle. 
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VII.  CONCLUDING  REMARKS 


The  following  conclusions  have  been  reached  as  a  result  of  this  study: 

a.  A  time -dependent  technique  with  the  MacCormack  finite  difference 
scheme  provides  stable  integration  for  both  one-  and  two-phase 
nozzle  flow  equations. 

b.  The  utilization  of  the  BFC  system  enhances  the  capability  of  the 
program  to  the  solution  of  flow  inside  nozzles  with  complex 
geometry. 

c.  Imbedded  shock  can  occur  in  the  region  downstream  of  the  nozzle 
geometric  throat  for  the  flow  inside  nozzle  with  small  throat 
radius  of  curvature. 

d.  The  small-sized  particles  act  more  effectively  to  slow  down  the 
gas -phase  expansion  than  that  of  large- sized  particles  for  the 
same  particle  mass  fraction. 

e.  For  a  two-phase  flow  with  high  particle  loading  ratio,  the  gas- 
phase  can  become  subsonic  at  the  geometric  throat. 

f.  The  computed  one-  and  two-phase  results  are  important  for 
nozzle  wall  heat  transfer  and  ablation  study. 

g.  In  general,  the  assumption  of  constant  fractional  lag  is  not  justi¬ 
fied  for  a  two-phase  transonic  flow.  The  prediction  of  the  gas- 
particle  flow  field  requires  that  the  proper  momentum  and  energy 
exchange  between  the  gas  and  particles,  such  as  the  fully  coupled 
solution  presented  in  this  study,  be  taken  into  account. 
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